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Q 
U 

' We study classical and quantum maps on the torus phase space, in the presence 

of noise. We focus on the spectral properties of the noisy evolution operator, and 
prove that for any amount of noise, the quantum spectrum converges to the classical 
one in the semiclassical limit. The small-noise behaviour of the classical spectrum 
highly depends on the dynamics generated by the map. For a chaotic dynamics, the 
I outer spectrum consists in isolated eigenvalues ("resonances") inside the unit circle, 

leading to an exponential damping of correlations. On the opposite, in the case of 
a regular map, part of the spectrum accumulates along a one-dimensional "string" 
, connecting the origin with unity, yielding a diffusive behaviour. We finally study 

2 I the non-commutativity between the semiclassical and small-noise limits, and illus- 

trate this phenomenon by computing (analytically and numerically) the classical and 
quantum spectra for some maps. 



Introduction 



> 

X 

■ Numerical studies of chaotic dynamical systems inevitably face the problem of rounding er- 
rors due to finite computer precision. Indeed, the unstability of the dynamics would require 
an infinite precision of the initial position, if one wants to compute the long-time evolu- 
tion. Besides, any deterministic model describing the evolution of some physical system is 
intrinsically an approximation, which neglects unknown but presumably small interactions 
with the "environment". A way to take this interaction into account is to introduce some 
randomness in the deterministic equations, for instance through a term of Langevin type. 
Such a term induces some diffusion, that is, some coarse-graining or smoothing of the phase 
space density. If the deterministic part of the dynamics is unstable, it has the opposite ef- 
fect to transform long- wave-length fluctuations of the density into short-wave-length ones, 
so the interplay between both components (deterministic vs. random) of the dynamics is 
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a priori not obvious. Rigorous results on the behaviour of noisy chaotic systems have been 
obtained recently j23IEllH|: they show that a usual property of deterministic chaos, namely 
the exponential mixing, still holds in presence of noise. Actually, introducing some noise 
provides a way to compute the exponential decay rate (and the subsequent "resonances") 
in a controlled fashion P EEl ED ESI Ei ■ 

In quantum mechanics, a (small) system is never perfectly isolated either, and is subject 
to inevitable interaction with the environment, responsible for decoherence effects. The 
study of decoherence has received a large attention recently, due to the (mostly theoretical) 
interest in quantum computation |T2| and precise experiments measuring decoherence. In 
some simple cases, one can study the dynamics of the full system (small plus environment), 
and obtain an effective dynamics of the small one [H]. Under suitable assumptions, this 
effective dynamics results in a memoryless (super) operator acting on the quantum density 
of the small system. 

A possible way to modelize this operator is by quantizing the "Langevin term" used 
in the classical framework, that is introduce a random noise in the quantum evolution 
equation (for either continuous or discrete time dynamics). Averaging over the noise yields, 
in the continuous-time framework, a non-Hermitian "Lindblad operator" [S^, while in the 
discrete-time case the evolution can be cast into the product of a unitary (deterministic) 
evolution operator with a "quantum coarse-graining" operator this product being called 
from now on the "coarse-grained evolution operator". Similar operators have also appeared 
in the theoretical study of spectral correlations for quantized maps [SS], and in models of 
dissipative quantum maps on the sphere |18| . 

Several recent studies have been devoted to noisy or coarse-grained quantum maps 
[UlEIllZl in the semiclassical limit, emphasizing the respective roles of regular vs. chaotic 
regions in the classical phase space. The aim was either to study the time evolution of 
an initial density or to compute the spectrum of the classical or quantum coarse- 
grained evolution operators |26( I31| . A different type of quantum dissipation operator, 
originating from the theory of superradiance, was composed with a unitary quantum map 
on the sphere, and studied in a series a papers by D. Braun [T3j. The author computed 
a Gutzwiller-like semiclassical trace formula for powers of the full (non-unitary) propa- 
gator, and showed their connection with traces of the corresponding classical dissipative 
propagator. 

In the present paper, we present some results about the spectrum of coarse-grained 
propagators, for maps defined on the 2-dimensional torus (Section HI . For the sake of 
simplicity, we restrict ourselves to either fully chaotic or fully regular maps. In the limit 
of vanishing noise, the spectrum of the classical coarse-grained evolution operator behaves 
differently in these two (extreme) cases: for a chaotic map, the spectrum has a finite gap 
between unity and the next largest eigenvalue, due to exponential decay of correlations 
|ifl|l8]. while in the regular case some eigenvalues come arbitrary close to unity. To illustrate 
these results, we study in detail some linear systems (either chaotic or integrable), for which 
the eigenvalues may be computed analytically. 

We then turn to the quantum version of these systems (Section El) • After recalling the 
setting of quantum maps on the torus, we define the quantum coarse-graining operator. 
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and then prove that for any classical map and a fixed finite noise, the spectrum of the 
quantum coarse-grained evolution operator converges to the spectrum of the classical one 
in the semiclassical limit (Theorem QJ. 

We finally study in Section |31 the non-commutativity of the semiclassical vs. vanishing- 
noise limits, using as examples the maps studied in the classical framework (Section HJ. As 
a byproduct, we show that one can obtain the resonances of a classical hyperbolic map from 
the spectrum of an associated quantum operator (the quantum coarse-grained propagator), 
provided the coarse-graining is set to decrease slowly enough in the semiclassical limit (we 
conjecture a sufficient condition for the speed of convergence). On the opposite, for an 
integrable map the same limit yields a spectrum densely filling one-dimensional curves 
("strings") in the unit disk, one of them containing unity as a limit point. 

1 Classical noisy evolution 

The classical dynamical systems we will study are defined on a 2-dimensional symplectic 
and Riemannian manifold, the torus = M^/Z^. The maps are smooth (C°°), invertible 
and leave the symplectic form dp A dq (and therefore the volume element d'^x = dq dp) 
invariant: they are canonical diffeomorphisms of T^. Such a map x = {q,p) ^ Mx 
naturally induces a Perron- Frobenius operator V = Vm acting on phase space densities 
p{x): [Vp]{x) = p{M~^x). 

1.1 Spectral properties of the classical propagator 

In this section we review the spectral properties of the Perron- Frobenius operator Vm, 
depending on the map M and on the functional space on which the operator acts. The 
results presented are not new, but allow us to fix some notations. 

1.1.1 Spectrum on ^^(T^) 

For any canonical diff'eomorphism M and any p > 1, the spectrum of Vm on the Banach 
space LP(T^) is a subset the unit circle and admits ^0(2^) = 1 as invariant density. In 
particular, V is unitary on L^(T^) and on its subspace LKT"^) of zero-mean densities. 

One can relate some dynamical properties of the map M with the (unitary) spectrum of 
V on Lq(T^) jl8j . For instance, if the dynamics M leaves invariant a nonconstant observable 
H{x) G C°(T^) (e.g. if M is the stroboscopic map of the Hamiltonian fiow generated by H), 
then all observables of the type p{x) = f{H{x)) (with / a smooth function) are invariant 
as well, so that the eigenvalue 1 of P is infinitely degenerate. The full spectrum of V will 
be explicitly given for some integrable maps in Section 12.21 

On the opposite, the map M is ergodic iff V has no invariant density in Lq (i.e. if 1 
is a simple eigenvalue in L^). Stronger chaotic properties may be defined in terms of the 
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correlation function between two densities f^gELi^: 

Cfg{t) [ dxf{x)g{M-'x) = {f,V'g) (1) 

(the time parameter t will always take integer values). The map M is said to be mixing 
iff for any f,g, the correlation function behaves for large times as: 

^/g(^) '"^ — ^ / dxf{x) i dxg{x). (2) 

A slightly weaker property (weak mixing) is equivalent with the fact that Vm has no 
eigenvalue in Lq(T^). 

1.1.2 Exponential mixing 

For a certain class of maps (e.g. the Anosov maps defined below), the convergence of 
mixing is exponentially fast, provided the observables f,g are smooth enough. One speaks 
of exponential mixing with a decay rate 7 > if \Cfg{t) — J f dx J g dx\ < -ft'/^e"'^!*! for a 
certain constant Kjg. 

This behaviour can sometimes be explained through the spectral analysis of the Perron- 
Frobenius operator V acting on an ad hoc functional space B, with B L"^ (in general 
B contains some distributions). One proves that the operator V oiv B is quasi- compact: 
its spectrum consists of finitely many isolated eigenvalues {Aj} (called Ruelle-Pollicott 
resonances) situated in an annulus {r < Aj < Aq = 1} for some r > 0, plus possible 
essential spectrum inside the disk of radius r. In that case, assuming for simplicity that 
each Aj is a simple eigenvalue with spectral projector Ilj, the spectral decomposition of V 
on B leads to the following asymptotic expansion for Cfg(t) in the limit t ^ 00: 

Cfg{t) = Y,Xl{f,U,g) + 0{r'). (3) 

i 

In the above expression, the brackets do not refer to a Hermitian structure, but represent 
the integrals / f{x)[Ilig]{x)dx. The first resonance Aq = 1 corresponds to the unique 
invariant density po, and the decay rate (obviously independent of the observables f,g) is 
given by 7 = — log |Ai|, resp. by — logr if there is no resonance other than unity. 

This type of spectrum was first put in evidence in the case of uniformly hyperbolic 
maps by using Markov partitions to translate the dynamics on the phase space into a 
simple symbolic dynamics (namely a subshift of finite type) jlOI. It was later extended to 
more general systems, including non-uniformly hyperbolic ones In the next sections, 
we will introduce the Anosov diffeomorphisms on the torus, which often serve as "model" 
for deterministic chaos. 

1.1.3 Anosov difFeomorphism on the two-dimensional torus 

In this section we recall the definition and some properties of an Anosov difFeomorphism M 
on the torus |2l|. The Anosov property means that at each point x G T^, the tangent 
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space TxT'^ splits into T^T'^ = © i?", where Ex are the local stable and unstable 
subspaces. The tangent map d^M sends Ex^^ to E^^^^,^^ and there exist constants A > 0, 
< As < 1 < A„ (independent on x) s.t. 

Vt e N, ||(4M*)|s.|| < AXl and || (4M-*)|s. || < AX-'. (4) 

These inequalities describe the uniform hyperbolicity of M on T^. The splitting of T^-T^ into 
Ex^'^ has in general regularity C^+° for some 1 > a > 0, meaning that it is differentiable 
and its derivatives are Holder-continuous with coefficient a. 

This uniform hyperbolicity implies that M is ergodic and exponentially mixing w.r.to 
the Lebesgue measure (see next section). 

Simple examples of Anosov diffeomorphisms are provided by the linear hyperbolic au- 
tomorphisms of T^, defined by matrices A G S'L(2,Z) with |try4| > 2. These maps 
are sometimes referred to as "generalized cat maps", in reference to Arnold's "cat map" 
/ 2 1 \ 

^Arnold = ( j [3|- Wc will study in some detail these linear maps and their quantiza- 

tions, and obtain a rather explicit description of the associated coarse-grained propagators. 
One can perturb the linear hyperbolic automorphism A with the stroboscopic map (p\j gen- 
erated by some Hamiltonian H{x) on the torus: M =^ (p\j o A. If the Hamiltonian H is 
"small enough", the perturbed map M remains Anosov, and topologically conjugated with 
the linear map A. 



1.1.4 Ruelle resonances for Anosov diffeomorphisms 

As announced above, Anosov diffeomorphisms on are exponentially mixing, and their 
correlation functions satisfy expansions of the type Although the original proofs made 
use of Markov partitions we will rather describe a more recent approach due to Blank, 
Keller and Liverani j^, which has the advantage to provide spectral results for the coarse- 
grained propagator as well. 

The strategy of |H1 is to define a Banach space B of densities on adapted to the map 
M: the space B contains distributions which are smooth along the unstable direction of 
M but possibly singular (dual of smooth) along the stable direction. In particular, these 
functions are too singular to be in L^, which allows a non-unitary spectrum of P acting on 
B. The space B is defined in terms of an arbitrary parameter < /3 < 1, and satisfies the 
continuous one-to-one embeddings: C^(T^) ^ B ^ C^(T^)*. 

The space B depends on the map M, and also on the direction of time: the space 
Bm-i adapted to the map is different from Bm- Although the map M on the torus is 
invertible, the dynamics of the operator V on the space B is irreversible: its spectrum is 
qualitatively very different from that of . 

The authors indeed show that V is quasicompact in B, with essential spectral radius 
Tess bounded above by a = max(A~^,Af) (see Eq. (jH) for the definition of A„/s). For 

any 1 > r > a, the spectrum of V in the ring Rr '= {r < |A| < 1} consists in isolated 
eigenvalues, the Ruelle-Pollicott resonances Therefore, a spectral expansion similar 
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with Q holds for any pair of observables f,gEB, which includes in particular observables 
in C^(T^) (the expansion might be slightly more complicated than Q due to possible finite 
degeneracies of the resonances). 

A possible strategy to extend the results of [H] to maps and observables of regularity C'' 
was discussed in the recent preprint [S]. Under stronger smoothness assumptions, namely 
for real-analytic Anosov maps in two dimensions, H. Rugh |41| constructed a transfer 
operator acting on observables real-analytic along the unstable direction, and "dual of 
analytic" along the stable one; he showed that this operator is compact, which means that 
the essential spectral radius vanishes in that case. 

1.2 Coarse-grained classical propagator 

In this section we precisely define the operator representing the effect of "noise" on the 
deterministic evolution of M. This operator is of diffusion type, it realizes a coarse-graining 
of the densities. 

1.2.1 Classical difFusion operator 

We consider a smooth probability density K{x) on with compact support. For sim- 
plicity, we also assume that K{x) = K{—x). From there, to any e > corresponds a 
probability density on the torus, defined as K^{x) = ^ Ylnei? ^ (2±ii). Due to the com- 
pact support of we see that for small enough e and any x G this sum has at most 
one nonvanishing term (for x close to the origin). We define the coarse-graining operator 
on L^(T^) as the following convolution: 

V/gL2(t2), [DJ]{y)= [ dxK,{y-x)f{x). (5) 

is a self-adjoint compact operator on L^(T^), with discrete spectrum accumulating at 
the origin. 

We define the Fourier transform on the plane as K{^) = dx K {x)e^™^^^ , with the 
wedge product given by ^ A x = ^ip — ^2<i- From the assumptions on the density the 
function K{^) is real, even and smooth. It takes its maximum at ^ = (where it behaves 
as K{^) = 1 — Q{C) + C'(l^l^), with Q{.) a positive definite quadratic form), and decreases 
fast for |,^| 00. 

The plane waves (or Fourier modes) on are accordingly defined as Pk{x) =^ exp {2mx A 
for k They obviously form an orthonormal eigenbasis of the coarse-graining operator: 

Ve > 0, VA; G Z^, D,pk = k{ek)pk. (6) 

The fast decrease of the eigenvalues as — > 00 implies that the operator is not only 
compact, but also trace-class. It kills the small-wave-length modes, effectively truncating 
the Fourier decomposition of p(x) at a cutoff \k\ ~ e^^. For some instances, we will actually 
replace the smooth function K{^) by a sharp cutoff 6(1 — |,^|) (with O the Heaviside step 
function). 
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1.2.2 Classical coarse-grained propagator 

The noisy dynamics associated with the map M is represented by the product of the 
deterministic evolution Vm with the diffusion operator: 

VM,e = D,oVM. (7) 

It may be more natural to define the noisy propagator as the more 'symmetric' D^oVm°D^, 
but both definitions lead to the same spectral structure. VM,e describes a Markov process 
defined by the deterministic map M followed by a random jump on a scale e. 

We now give some general properties of this operator, independent of the particular map 
M. Like the regularizing operator D^, 'Pm,^ maps distributions into smooth functions, and 
is compact and trace-class on any functional space containing C°°(T^) as a dense subspace, 
with a spectrum independent on the space. Its eigenvalues {A^_e}/i>o are inside the unit disk 
(only Ao,e = 1 is exactly on the unit circle), they are of finite multiplicity and admit the 
origin as only accumulation point. The eigenvalue Ao,e is simple, with unique eigenfunction 
Po- T^e maps a real density to a real density, therefore its spectrum is symmetric with 
respect to the real axis. 

In the next section we investigate in more detail the behaviour of these eigenvalues in 
the limit of small noise, stressing the difference between chaotic vs. regular maps. 



2 Spectral properties of classical coarse-grained propa- 
gator 

We describe more precisely the spectrum of Pm,^, in the limit of small noise, and for 
different classes of maps M. We start with the most chaotic maps, namely the Anosov 
diffeomorphisms defined in Section 11.1.31 the exponential mixing of which was described 
in section IT. 1.41 In a second part, we will then turn to the opposite case of "regular" maps 
on T\ 



2.1 Anosov difFeomorphism 



We use the same notations as in Section 11.1.41 for M an Anosov map. It was proven in 
[H] that the spectrum of VM,e outside some neighbourhood of the origin converges to the 
resonance spectrum of Vm on the Banach space B. 

More precisely, for M a smooth Anosov canonical diffeomorphism on with foliations 
of Holder regularity C^"*"" (0 < a < 1), one considers the "associated" Banach space B = Bm 
defined in terms of a coefficient P < a, such that the essential radius of P on i3 has for upper 
bound a = max(A~\ Af ). The authors then construct a "weak" norm || . on C{B), such 
that llPe — PIL — ^ as e — s> 0. As a consequence, for any 1 > r > a, any A in the annulus 

n 



{r < |A| < 1} and any 6 > small enough, the spectral projector H^j^^^) of V, 
of V) in the disk B{X,6) = {z : \z — X\ < 6} satisfy 



B{X,S) 



n 



(^) 

B{X,S) 



n 



B{X,S) 



(resp. 
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Both projectors therefore have the same rank for e small enough, this rank being zero if 
the disk contains no resonance Aj. 

This proves that the spectrum of in the annulus Rr converges (with multiplicity) 
to the set of resonances {Aj} as e — > 0, and the eigenmodes of Ve weakly converge to 
corresponding eigenmodes of V (since the latter are genuine distributions, the convergence 
can only hold in a weak sense). 

Remarks: 

• These results also apply to the 'symmetric' coarse-graining D^oV o D^. 

• It is reasonable to conjecture that these results hold as well if the coarse- graining 
kernel K{x) is not compactly supported on M^, but decreases sufficiently fast, for 

instance if one takes the Gaussian G{x) =^ e^'^'^'l^, G{^) = e^''!^!^, as was done in 
[H [2^1. As a broader generalization, we will sometimes consider a coarse-graining 
defined by a sharp cut-off in Fourier space, K{^) = 6(1 — |^|), similar to the method 
used in IHl]; a finite-rank coarse-graining was also used in for 1-dimensional noisy 
maps. 

• If the map M and the kernel K{x) are real-analytic, we conjecture that the eigen- 
values of on any ring Rr with r > converge to the resonances V in Rr (cf. the 
remark at the end of Section fl .1.411 . 

The (discrete) spectrum of Ve is the same on any space S admitting C°° as dense 
subspace, in particular on L^, and this is the space we will consider from now on (more 
precisely, its subspace Lq). This spectrum drastically differs from the absolutely continuous 
unitary spectrum of the "pure" propagator V on that space (cf. Section ITTT|l . The unitary 
spectrum is thus unstable upon the coarse-graining D^: when switching on the noise, the 
spectral radius of Ve on Lq suddenly collapses from 1 to |Ai| < 1. As we will see in the 
next sections, this collapse is characteristic of chaotic maps. We first describe it explicitly 
for the case of hyperbolic linear automorphisms. 

2.1.1 Example of Anosov maps: the hyperbolic linear automorphisms 

In this section we review [THl Ell the spectral analysis of the (pure vs. noisy) propagator 
when the map A is a hyperbolic linear automorphism of (cf. Section H.l. 3 jl . The unitary 
operator Va on Lq acts very simply on the basis of Fourier modes pk, k E Zl = Z"^ \0, 
namely as a permutation: 

[PAPk] (x) = pk{A'^x) = pAk{x). (8) 

This evolution induces orbits on the Fourier lattice Z^, which we will denote by O{ko) = 
{A^ko, t E 1,}. Due to the hyperbolicity of A, each orbit is infinite, so that the modes 
{Pyi*A:o) t E Z} span an infinite-dimensional invariant subspace of Lq, which we call SpanOiko). 
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The spectral measure of Va associated with this subspace is of Lebesgue type (as usual 
when the operator acts as a shift [39|). The number of distinct orbits being infinite, the 
spectrum of Va on Lq is Lebesgue with infinite multiplicity. 

Now we consider the coarse-grained propagator Ve,A = o Va- Since the pk are 
eigenfunctions of (cf. Eq. (jH])), Ve,A will also act as a permutation inside each orbit 
(9(/co), but now at each step the mode pk is multiplied by K{ek). Therefore, the operator 
Ve,A restricted to the invariant subspace SpanO^ko) can be represented as follows (the 
notation (., .) denotes the scalar product on Lq): 

Ve,A\oiko) = ^K{eA^^^ko) {pA^ko, ■ ) Pa'+i/cq- (9) 

Since |v4*/co| — ^ oo for t ^ ±00, the factors K^eA^^^ko) vanish in both limits. As a result, 
the spectrum of ^£,^10(^0) reduces to the single point {0} (which is not an eigenvalue, but 
essential spectrum) [39J. By taking all orbits into account, the spectrum of V^^a on Lq also 
reduces to {0}: for linear hyperbolic maps, the collapse of the Lebesgue unitary spectrum 
through coarse-graining is 'maximal'. 

2.2 Coarse-graining of regular dynamics 

In the previous sections we have considered classical propagators of Anosov diffeomor- 
phisms on T^. We now describe the opposite case of an integrable map on T^. The notion 
of integrability for a map is not so clear as for a Hamiltonian flow. In view of the examples 
below, the definition should include the stroboscopic map = ip\j of the flow generated 
by an autonomous Hamiltonian H[x) on T^, but it should also encompass elliptic and 
parabolic automorphisms, as well as rational translations. A required property is that the 
phase space splits into a union of invariant 1-dimensional (not necessarily connected) 
closed submanifolds, with possible "critical energies". Note that this condition excludes the 
(un)stable manifolds of an Anosov map, which are open. A more or less equivalent condi- 
tion for integrability is that the map leaves invariant a nowhere constant smooth function 
H{x), the level curves of which provide the above submanifolds. As a result, any density 
p{x) = f{H{x)) is invariant through Vm as well, so that Vm has an infinite-dimensional 
eigenspace of invariant densities on Lq (which we will call Vinv)- 

The rest of the spectrum of Vm can be of various types, as we will see (it can be pure 
point or absolutely continuous, be a mixture, etc.). For this reason, a general statement 
concerning the coarse-grained spectrum of these maps cannot be very precise. 

In the following sections we consider some simple examples of integrable maps for which 
(part of) the spectrum can be analyzed in detail. We then discuss (mostly by hand- waving) 
the general case. 

2.2.1 Translations on the torus 

The simplest nontrivial maps on are the translations x T^x = x + v mod T^, which 
do not derive from a Hamiltonian on the torus. A translation is either ergodic (yet non- 
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mixing) or integrable (see below). 

The spectrum of the corresponding Perron-Frobenius operator Vv = Vt^ on Lq is easy 
to describe (TH]: Vv admits the Fourier modes pk as eigenstates, with eigenvalues e^''^'^'^^. 
The spectrum of Vv is therefore pure point, with possible degeneracies. The spectrum of 
the coarse-grained propagator Pe^„ = D^oVv on Lq is also easy to describe: each Fourier 
mode Pk is an eigenfunction with eigenvalue K[ek)G^™^^^ inside the unit disk. From the 
small-e expansion K{ek) ~ 1 — e^Q{k), the eigenvalues corresponding to long wavelengths 
(|fc| <C e~^) are close to the unit circle for small e, while the short wavelength eigenvalues 
(|/c| ^ e^^) accumulate near the origin. We now describe how the global aspect of the 
spectrum qualitatively depends on the translation vector v. 

• Tv is ergodic iff the coefficients fi, V2 as well as their ratio ^ are irrational: in that 
case, the equation A G Z has no solution for k ^ The spectrum of Vv forms a 
dense subgroup of the unit circle, all eigenvalues being simple. In the limit e 0+, 
the spectrum of V^^v becomes "dense" in the unit disk (a similar quantum spectrum 
is plotted on Fig. |31 right). 

• If one of the coefficients f i, f2, ^ is rational, Tv leaves invariant a family of parallel 
1-dimensional "affine" submanifolds, and is therefore integrable according to our def- 
inition. For instance, if we take ^^i = ^ (with r, s coprime integers) and V2 irrational, 
then for any go ^ [0, 1], the union of vertical lines |Jz*=o ~ ^0 + ^} is invariant (if 
s > 1, this set is non-connected). The spectrum oiVv is still dense on the circle, but 
all eigenvalues are now infinitely degenerate: for any /cq, the modes A; = /cq + (0, js), 
j G Z share the eigenvalue _ 'jj^g eigenvalues of V^^v are at most finitely- 
degenerate: to each phase e^''^''°^^ corresponds a "string" of eigenvalues of decreasing 
moduli, the largest one being at a distance ~ from the unit circle. As in the 
previous case, the spectrum densely fills the unit disk as e — 0. 

• If both fi,f2 are rational of the form ^ with gcd(ri, r2, s) = 1, each point of 

is periodic with period s, the map is integrable. The only eigenvalues of Vv are of 
the phases e^''"-'/^, all being infinitely degenerate. The eigenvalues of Ve,v are at most 
finitely degenerate, they are grouped into s strings of phases e^™^^^, j = 0, . . . , s — 1. 
For small e, the eigenvalues become dense along these strings, the largest eigenvalue 
at a distance ~ from the unit circle (this spectrum is similar with the quantum 
one plotted on Fig. HJ left). 

Comparing the first and second case, we see that the spectrum ofV^ cannot unambiguously 
differentiate an ergodic from an integrable map: both may have eigenvalues close to the 
unit circle. On the other hand, the second and third cases both correspond to integrable 
maps. Yet, their small-noise spectra look quite different from one another. 

2.2.2 Non-hyperbolic linear automorphisms of the torus 

Another class of non-mixing linear maps on is provided by the non-hyperbolic linear 
automorphisms. These automorphisms split into two classes (for notations, we refer to 
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Section EmU): 



the elliptic transformations (|trA| < 2), like for instance the 7r/2-rotation given by 
the matrix J = ( ^ J" ] . As opposed to the hyperbolic case, each Fourier orbit 



(9(fco) = {J-'ko, j = 0, ...,3} is finite of period 4, and SpanO^ko) splits into 4 
eigenspaces associated with the eigenvalues {i^, / = 0, ... ,3}. Switching on coarse- 
graining, the eigenvalues of Vej read i'Afcg, with A^^ = \J\K{tkQ)K{tJkQ) \ . The 
spectrum of 'Pe,j on Lq therefore consists in four strings along the four half-axes, which 
become dense in the limit e — >^ (see fig. [T]for the analogous quantum spectrum). 
Similarly, an elliptic transformation of trace iiA = 1 will satisfy = 1, therefore 
the spectrum of P^.A forms 6 strings. An elliptic transformation of trace tiA = — 1 
will lead to 3 strings. 

the parabolic transformations (or parabolic shears), given by matrices of the type 

5'=(q ^ j , s E Z^. The dynamics reads {q,p) ^ {q + sp,p), so any periodic 

function of the momentum p is a conserved quantity. The Fourier vector ko = {ki, ^2) 
generates the orbit C(fco) = {ko + {jsk2,0), j G Z}: if ^2 = 0, the mode pk^ is 
invariant and the orbit is a singleton; on the opposite, if k2 7^ 0, Oi^ko) is infinite, and 
leads to Lebesgue spectrum for Vs\o(ko)- The full spectrum ofVs on Lq is therefore the 
union of the infinitely degenerate eigenvalue 1 with a countable Lebesgue spectrum 
[T8] . In the case ^2 7^ 0, the noisy propagator Ve,s acts on SpanO{ko) as in Eq. 
upon replacing A by S. Since the wavevectors \S^ko\ diverge in both limits t ±00, 
the associated spectrum reduces to the singleton {0}. On the opposite, each mode 
P(ifi) is eigenstate of V^^s with eigenvalue K (e(/, 0)): these modes form a string along 
the real axis, which becomes dense as e — 0. 



2.2.3 Nonlinear shear 

If we perturb the linear parabolic shear S (cf. last section) with the stroboscopic map 
of the flow generated by a Hamiltonian of the form F{p), we obtain a nonlinear shear 
{q,p) ^ {q + sp + F'{p),p). We assume that s > 0, and that the Hamiltonian F satisfles 
s + F"{p) > everywhere. This map still conserves momentum, so that any density p{p) 
is invariant. Among these, the Fourier modes p(/ o) are eigenstates of "P^, with eigenvalues 
i^'(e(/,0)): they form the same string of real eigenvalues as for the linear shear. 

For any n E Z^,V and leave invariant the subspace Vn = Span |e^'™'^p(p), p G L^(T^) }, 
and their spectra on this subspace can be partially described, at least if one takes F' small 
enough and replaces the kernel K{^) by the sharp cutoff 6(1 — |^i|)6(l — |^2|) (such that 
projects to a flnite-dimensional subspace). The spectrum of P on 14 is absolutely contin- 
uous, while that of Ve is contained in a small neighbourhood of the origin, uniformly with 
respect to n and e (see Appendix lA. II for details). The spectrum should be qualitatively 
the same if i^(^) is of fast decrease at infinity. 
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2.2.4 Stroboscopic map of the Harper Hamiltonian flow 

As a last example of integrable map, we consider the Harper Hamiltonian on the torus 

H{x) = cos(g) + cos(p) = ^(p(i,o) + P{-i,o) + P(o,i) + P(o,-i)), 

and take its stroboscopic map M = ip\j. The invariant densities are of the form p{x) = 
f{H{x)). As opposed to what happened for the parabolic shear, Vinv is not invariant 
through the diffusion operator D^, so the spectrum of "P^ is more complicated to analyze. 
Still, if we take a coarse-graining kernel satisfying K{C,) = A; (|^i| + |^2|) along the axes of 
slopes 0, ±1/2, ±1, ±2, oo) then the invariant functions H{x), {H{x))'^ — 1 and {H{x))^ — 
are eigenstates of and therefore of Ve, with respective eigenvalues fc(e), k{2e) 
and k{3e) (if the kernel only satisfies K{{0,()) = -^((C)O)) together with even parity, then 
H{x) will be eigenstate of with eigenvalue i^((e, 0))). These eigenvalues are real and 
approach unity as e ^ 0. Numerically, they are the three first elements (resp. the first 
element) of a string of real eigenvalues connecting unity to the origin, the further elements 
of the string being mixtures of invariant and non-invariant densities (see Fig. EJ right for a 
similar quantum spectrum). 

2.2.5 General behaviour for an integrable map 

After these examples, we want to describe the spectrum of VM,e for M an integrable map, 
say the stroboscopic map of an autonomous Hamiltonian H. As we already explained, the 
space Vinv of invariant densities is infinite-dimensional, containing all densities f{H{x)). 
When applying the coarse- graining D^, one generally mixes these invariant densities with 
non-invariant ones (as opposed to what happens for the linear maps described above). 
Using ideas of degenerate perturbation theory, we conjecture the following behaviour, which 
is supported by numerical investigations |3T| l26]. 

The invariant subspace Vinv contains modes piw with long- wavelength fluctuations (that 
is, densities for which the Fourier spectrum is concentrated in a finite region near k = 0). 
For e small, these modes are hardly modified by the coarse-graining operator: D^piw — pivj = 
0{e^). This suggests that has an eigenstate of the type pi^ + C'(e^), of eigenvalue 
1 — O(e^), this eigenstate being close to an invariant state of V. On the opposite, Vinv 
also contains highly fluctuating modes, which will be very damped by D^, and should lead 
to eigenstates of Ve close to the origin. As a whole, the hybridization of invariant modes 
with noninvariant ones leads to a string of eigenvalues connecting unity to the origin, 
this string being symmetrical w.r.to the real axis. Notice that the projected operator 
^invVJlinv = ^invDJlinv = (whcrc Hinv IS the orthogoual projector on the space Vinv) is 
self-adjoint, therefore its spectrum can be estimated through the min-max method: for 
e — * it consists in a dense string of eigenvalues between unity and zero. I believe that the 
largest eigenstates (-values) of this projected operator are close to eigenstates (-values) of 
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Conclusion: qualitatively different spectra 

We have exhibited a qualitative difference of the coarse-grained Perron-Frobenius spectra 
between, on one side, the chaotic maps (Anosov), on the other side, non-mixing maps, 
including ergodic translations and integrable maps. 

In the latter case, the spectrum of Ve on the subspace Ll{T'^) comes close to the unit 
circle for small e, either along 1-dimensional strings "connecting" the origin to infinitely 
degenerate eigenvalues of V (among which unity), or by densely filling the unit disk (for 
ergodic translations). A common feature is that any annulus {R < \\\ < 1} (or any open 
neighbourhood of unity) contains more and more eigenvalues of Ve in the small-noise limit 
e ^ 0. 

On the opposite, for an Anosov map the spectrum of Ve on Lq is contained inside 
a disk of radius R < 1, uniformly for small enough e. This spectrum consists in a finite 
number (possibly zero) of finitely degenerate eigenvalues (asymptotically close to the Ruelle 
resonances) inside an annulus {r < |A^e| < R}, the remaining eigenvalues having moduli 
smaller than r. 

The numerical results of l^llEEl go beyond this statement: the authors consider systems 
with mixed dynamics (not to be mistaken with the "mixing" property of Anosov maps), that 
is systems for which the phase space splits into "regular islands" embedded in a "chaotic 
sea". They study the spectrum of the coarse-grained Perron-Frobenius operator (the coarse- 
graining is taken as a cut-off in Fourier space), and show the presence of eigenvalues close 
to the unit circle, the eigendensities of which are supported on the regular islands; on the 
other hand, they also find some eigenvalues inside the unit circle, which are associated 
with the chaotic part of phase space, and therefore interpreted as (generalized) Ruelle 
resonances. Yet, the nature of the "chaotic sea" in such systems is not well understood at 
the mathematical level, so that a rigorous spectral analysis of the propagator seems a quite 
distant goal. 

3 Quantum coarse-grained evolution 

After studying the effect of noise on classical propagators, we turn to noisy quantum maps, 
obtained by quantizing the classical ones. Although we will restrict ourselves to maps on 
the 2-torus, the main result of this section (TheoremHJ can be straightforwardly generalized 
to quantum maps on any compact phase space, provided one adapts the definition of the 
coarse-graining operator (see the discussion in Section ri2.3|) . We start by recalling the 
setting of quantized maps on the 2-torus. 

3.1 Quantum propagator on the torus 
3.1.1 Quantum Hilbert space and observables 

We shortly review the construction of quantum mechanics on T^, in order to fix notations. 
For any value of /i > 0, the Weyl-Heisenberg group associates to each vector v = (fi,f2) G 
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the "quantum translation" Ty = exp |-^(t>2g — tiip)} which acts unitarily on L^(]R). 

These translations satisfy the group relations T^T^i = e 

The "torus wavefunctions" are then defined as distributions \ip) G iS'(R) satisfying 
the periodicity conditions T(o,i)|'?/') = T^i,o)\'^) = 1^)- Due to the group relations, such 
distributions exist iff h satisfies the condition {27ih)~^ = N E N (such a value of h will be 
called admissible). In that case they form a space of dimension A^, which will be noted Hn 
[2ni- A basis of this space is provided by the "Dirac combs" {\(lj)N}j=o Ar_i defined as: 

{Q\Qj)N = ^S{q-qj-u), with gj = (10) 



By construction, \qj) = |gj+7v), so the index j must be understood modulo A^. For prac- 
tical reasons, we will choose the following representatives of integers modulo A^: = 
{-f + 1, . . . , f } for iV even, = {-^ + 1, • • • , ^} for N odd. 

The quantum translation acts inside 7i at iff f is on the square lattice of side 1/iV, 
that is f = (^, ^) with G Z. Besides, each translation with integer coefficients acts on 

Hn as a multiple of the identity. As a result, the set < Ty/N, forms a basis of the 

space of linear operators on TCn, denoted by C{TCn)- This basis can be used to define the 
Weyl quantization of smooth observables / G C°°(T^). From the Fourier decomposition 
/ = J2k&'2 f{k)pk, the Weyl quantization of the observable / is defined as 

/ = Op;v(/) = $^/»T,/^= Yl ^^/^ [Y.i-l)''''^^'^'''^ hk + Nu) \ . (11) 

The "converse" of Weyl quantization, that is the Weyl symbol (or Wigner function) 
PFg(x) of an operator B G C(Hn) is also easily defined through its Fourier transform: 

for each k G Z^, its Fourier coefficient Wj^{k) is given by Wj^{k) = -^tr {tIi^B^. These 

Fourier coefficients are A^-periodic (up to a sign), so that the "function" Wj^{x) is a periodic 
combination of Dirac peaks on the lattice of side 1/2A^ jS^. Alternatively, a polynomial 
Weyl symbol was defined in [THj as the finite sum 



2i7rxAk 



As opposed to the former symbol, the polynomial symbol depends on the specific choice for 
the representative Z^r (the choice we made, with maximum symmetry around the origin, 
seems more natural in this respect). The polynomial symbol map together with Op at 

yield an isometric isomorphism between the subspace Xjy *== Span {p^, k G Z^} of L^(T^) 
and the space of observables on Hn equipped with the Hilbert-Schmidt scalar product 
{13, C) = -^tr(i?^C*). Since the Hilbert-Schmidt norm differs from the usual operator norm 
on C{Hn), we will denote by the space of observables on Hn (that is, A^ x A^ matrices) 
equipped with the Hilbert-Schmidt norm. 
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3.1.2 Quantization of canonical maps 

We briefly explain how one quantizes a canonical map M on T^. The aim is to deflne for 
each G N a unitary operator on Hn (which we will denote by Mtv or simply M) which 
satisfies prescribed semiclassical properties |12]- These properties do not unambiguously 
define the sequence of unitary matrices, so a choice has to be made (in a Toeplitz 
quantization is proposed for symplectic maps on Kahler manifolds). We present here 
another quantization prescription, which uses the following decomposition of any canonical 
map M on 

M = AoT^o(p]^. (12) 

In this formula, the linear automorphism A E SL{2, Z) and the translation Ty are uniquely 
defined. On the opposite, the last factor correspond to the stroboscopic map of the fiow of 
a time-dependent periodic Hamiltonian H G C°°(T^ x T^); this Hamiltonian is not unique, 
since two Hamiltonians Hi ^ H2 may lead to the same stroboscopic map = (p]j^ . From 
this decomposition, on can quantize M as follows |2^. First, one quantizes a la Weyl the 
Hamiltonian H{x,t) into the operator H{t) on Hn, then take for the quantization of ip]j 
the time-ordered exponential Te~^^o 'itH(t)/ft Second, one may quantize the translation 
Ty with the quantum translation T^(jv), where the vector f*^^) belongs to the 1/A^-lattice 
and is close to v, for instance take f^^^ = (^^'"iM^'"^]) Third, provided A satisfies 

the condition A = Id2 mod 2 or A = ax mod 2, the linear automorphism A is "naturally" 
quantized into a unitary matrix Aj^ [2^] (this condition may be relaxed if one generalizes 
the quantum Hilbert space Ti^ to arbitrary "Bloch angles" jl2j). Finally, the map M is 
quantized on Hn as the product: 

Mn = Ano f,(iv) o ^0 dtH(t)_ (13) 

This choice for the quantum map automatically satisfies the Egorov property [121 [221: for 
any classical observable p G C°°(T^), 

Mn Opn{p) - OpNip o M-i) 0. (14) 

This means that the quantization and finite-time evolution of densities commute in the 
semiclassical limit (the convergence holds for the operator norm in C{T-Cn))- 

In practice, the classical maps we consider are all defined as products of automorphisms, 
translations and Hamiltonian maps, so they admit a "natural" quantization. 

3.1.3 Propagator of quantum densities 

The quantum map M propagates quantum states {ip) E Hn- A density matrix p is an 
element of the space L% = Hn ^ '^at, and its evolution through the quantum map reads 

p MpM~^ = ad(M)p. The operator Vm '= ad(M) acting on the space is the 
quantum analogue of the Perron-Frobenius operator Vm acting on classical densities. Vm 
is unitary, with eigenvalues {e'^^^j"^'-'; i,j = 1,...,N}, where {e^^j} are the eigenvalues of 
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the matrix M. As a contrast with its classical analogue, the space of invariant densities 
through Vm is at least iV-dimensional, since it contains all the rank-1 projectors 
where are the eigenstates of M. 

The operator V acting on densities is called a superoperator, to contrast with an operator 
acting ouTCn- Being the adjoint action of a unitary matrix, it conserves the purity of the 
density, which means that a pure density p = |^)(^| is mapped onto a pure density 
Vp = \ip'){ip'\. We notice that V conserves the trace of the density, that is the quantum 
counterpart of the total probability of the density p. 

3.2 Quantum coarse-graining operator 

As we remarked above, the spectrum of V on Lj^ is qualitatively different from that of V 
on L^(T^): the former has at least N invariant eigenstates, while the latter may have only 
one if the map M is ergodic. In Section II.IL we explained how the spectrum of V was 
sensitive to the functional space on which V acts. This is not any more the case in the 
quantum framework, since the propagator is a finite-dimensional matrix. Still, we showed 
in Section[2lan alternative way to obtain the (nonunitary) resonance spectrum of an Anosov 
map, namely by introducing some noise and studying the spectrum of the noisy propagator 
Ve = D^oV on the space L^(T^). This procedure can also be carried out at the quantum 
level, by first defining a "quantum coarse-graining" or "quantum diffusion" (super)operator 
De, use it to construct a quantum coarse-grained propagator P^, then study the spectrum 
of the latter on L^. To connect the classical and quantum frameworks, we will consider 
the semiclassical limit N ^ oo. 

3.2.1 Definition 

We define the coarse-graining superoperator by analogy with the classical one (Eq. 
The latter can be expressed as follows: 

V/eL2(T2), {DJ){x)= [ dvK,{v)p{x-v)= I dv K,{v) {V.p) (x), (15) 

where is the Perron-Frobenius operator for the translation T^. Since Ty is quantized 
on Hn into the unitary matrix T^{n), a natural way to define a quantum coarse-graining 
superoperator would be through the integral 

dvK^{v) Vy(N). 

For convenience, we prefer a slightly different definition. The map v ^ v^^'^ is constant on 
squares of side 1/A^, so the above integral reduces to a finite sum over V G , with each 
operator Vv/n multiplied by the average of over the corresponding square. being a 
smooth function, this average is semiclassically close to the value Kf^iV/N), therefore for 
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large the above integral is well approximated by the sum: 

F{e,N) 



iV2 



K,iy/N) ad {Tv/n) . (16) 



JV 



The prefactor F(e, A^) is needed to guarantee that conserves the trace, that is DePo = Po- 
This factor is easily expressed in terms of the 2-dimensional theta function: 

^^(eiV,C) = 5^i^(eiV(z/ + C)), Cell (17) 

For Gaussian coarse-graining K{x) = G{x) = e"'"'"'!^, this function reduces to the product 
of two classical Jacobi 1-dimensional theta functions. In the limit eA^ ^ 1, the function 
9k is peaked around the point C = 0, due to the fast decrease of K{^). Now, one easily 
checks that F(e, N) = 9k{€N, 0)~^ converges to K{0) = 1 in the limit eA^ — > oo. 

3.2.2 Spectrum of the quantum coarse-graining operator 

The spectrum of -D^ on Lj^ is as easy to analyze as that of (see Eq. 0). Using the 
group relation 

ad(f,)T;, = e-^T^t,/, (18) 

we see that for any k G Z^, the quantum translation Tk/N (= OpN{pk)) is eigenstate of Z)^ 
with eigenvalue 

^^''^ - m 2^/^^^/^^" - ^k(6N,o) • ^^^^ 

V (i£^ 

For fixed k and eA^ oo, one has the asymptotics d!ff^ = K{ek) + 0{{eN)~°') for any 

power a > 0, so that the eigenvalue of associated with Tk/N converges to the eigenvalue 
oW^ associated with its symbol pk- The estimate is sharper in the Gaussian case: 

WkeZ%, d%^ = e-'^' + 0{e-'^), (20) 

so that the relative deviations between classical and quantum eigenvalues are exponentially 
small, uniformly on the modes k G {|A;i|, \k2\ < y(1 — 5)} (with S > fixed). For the 
Gaussian case, both classical and quantum spectra are strictly positive, which is not true 
in general. The spectrum of for a Gaussian noise is plotted in Fig. [T] (left). 

If we replace k{^) by the sharp cutoff 9(1 - |^|), then = I for \k\ < e-\ SJJJ = 
otherwise: the coarse-graining truncates the expansion (fTI| . keeping only short wavevectors 
(modulo A^). In other words, truncates the Fourier series of the polynomial Weyl symbol 
W-', keeping only the coefficients \k\ < e~^. 

A Fourier cutoff was also used as definition for the quantum coarse-graining in [HI], but 
there the cutoff was applied to the Husimi function of p instead of its polynomial Weyl 
symbol. 
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3.2.3 Probabilistic interpretation and generalizations of coarse-graining 



In a different physical framework (one-dimensional quantum spin chains), Prosen recently 
defined a similar quantum coarse-graining by truncating the densities on finite-dimensional 
spaces IHE]- In this case, the quantum densities are decomposed into sums of spin operators 
acting on finite sequences of spins (e.g. ai{x)ai{x + l) ■ ■ •cri(x+/)). The truncation consists 
in keeping only those operators for which I < e~^. 

The RHS of Eq. (fT6|) expresses the superoperator in the Kraus representation, i.e. 
as a sum J2j &d(-Ej), where the operators Ej on Hn satisfy the condition E^Ej = Mat. 
Kraus superoperators conserve the trace and the "complete positivity" of density matrices 



The RHS of Eqs. (|15I16|1 may be interpreted as a random global jump for the density, 
which jumps at a distance v with a probability oc K^{v). D^p and its quantum counterpart 
represent the average over all these random jumps. Since or do not depend on 
time, they therefore represent the classical and quantum versions of a memoryless Markov 
process. 

In a different scope, a superoperator similar to was used in |H3 to study the spectral 
correlations of the quantum map M in the semiclassical limit. Coarse-graining was there 
interpreted as an average over a set of quantum maps close to identity. The phase space can 
be an arbitrary (quantizable) symplectic manifold, and the classical and quantum averages 
are generated by a finite set of Hamiltonians {-f^j}j=i,...,/; is a smooth probability kernel 
in / variables with compact support of scale e around the origin. The classical and quantum 
coarse-graining operators are defined as: 



This scheme is more general than what we have done on the torus: one does not need any 
group action on the phase space, but only a sufficient number of Hamiltonians. We recover 
our previous definition on the torus if we take for "Hamiltonians" the multivalued functions 
Hi = p mod 1, H2 = —q mod 1 (these functions are not well-defined on T^, but their fiows 
are). The qualitative spectral features of the classical coarse-grained propagator ^ for 
small e should not depend on the selected family of Hamiltonians {Hj}, as long as the 
second-order operator X]j=i(^J?j)^ elliptic (in the above case on the torus, this is the 
Laplace operator ^ + which is indeed elliptic) fI7\ . 

In the framework of quantum mechanics on the 2-sphere, a different type of dissipative 
superoperator Vr = oad(M) was considered in [E], where D^- is a quantum dissipation 
operator obtained by integrating a quantum master equation during the time r. This 
dissipation operator first appeared in the study of superradiance in atomic physics (Tn]. As 
a main difference with our noisy propagator P^, the (unique) invariant eigenmode of Vr is 
different from Pq. This corresponds to the fact that the corresponding classical propagator 
Vr (obtained from Vr by taking the semiclassical limit) does not leave invariant the uniform 




(21) 
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density po, but rather a more singular measure, which may supported either on a discrete 
set of points, or even on a more comphcated "strange attractor". 

3.3 Quantum coarse-grained propagator 

We will now study the quantum analogue ofVe, that is the coarse-grained quantum prop- 
agator Ve = De° Vm- Similarly as the classical propagator, maps a Hermitian density 
to a Hermitian density; as a result, its spectrum is symmetric w.r.to the real axis. Like D^, 
Pe is a Kraus superoperator, and therefore realizes a quantum Markov process: the quan- 
tum density first evolves through the deterministic map ad(M), then performs a random 
quantum jump through ad(Ty/7v), with a probability oc K^{V/N). 

Linear combinations of quantum translations were considered in [Jj as models for de- 
coherence on the quantum torus. The authors studied the evolution of densities through 
an operator similar with P^, by computing the time evolution of the 'purity' tr(/3^). They 
took for M the quantized Baker's map, which is fully chaotic, yet discontinuous on T^. 
More recently, similar computations were performed for smooth nonlinear perturbations of 
cat maps |2T] . 

The author of uses Gutzwiller-type trace formulas to estimate the traces of powers 
of the dissipative quantum propagator tr('P"), in the regime of n fixed and h — > oo; 
from them he shows that each trace converges to the corresponding trace of the classical 
dissipative propagator tr(P") (this trace being given by a sum over periodic orbits as well). 

In the following I will not use any trace formula, but more basic semiclassical and 
operator-theoretic techniques to compare quantum and classical coarse-grained propaga- 
tors. Ve is indeed spectrally similar with its classical counterpart Ve- It conserves the trace: 
VePo = po, but in contrast with its noiseless version V, it has for unique invariant density 
Po, all other eigenvalues (in number A^^ — 1) being inside the unit disk. As a non-classical 
property, destroys purity: the image of a pure state p = is not a pure state. 

The following theorem, which is the central result of this paper, relates more precisely 
the spectra of and Ve in the semiclassical limit: it states the "semiclassical spectral 
stability" of coarse-grained propagators. 

Theorem 1. For any smooth map M on the torus and any fixed e > 0, the spectrum of 
the quantum coarse-grained propagator = De° Vm on L\ converges in the semiclassical 
limit N oo to the spectrum of the classical coarse-grained propagator Ve = D^o Vm on 
L^(T^). For any r > 0, the convergence is uniform in the annulus Rr = {r < \\\ < 1}. 

To compare the classical and quantum propagators, we use the isometry between the 
subspace 2^ of L^(T^) and L^, induced by the Weyl quantization Op^ and its inverse 

(see section 13.1.11 for notations). Ve is then isometric to crN{Ve) '= oVeO Op^ ° Hj^ 
(IIj^ projects orthogonally L^(T^) onto Xjy). We therefore need to compare the operators 
o^NiVe) and Ve on L^(T^). The crucial semiclassical estimate is the following lemma. 

Lemma 1. The finite-rank operators crjv('Pj converge to Ve in the operator norm on 
L2(T2), in the limit N ^ oo. 
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Proof: The key semiclassical ingredient is Egorov's property ()14|1 . From the norm 
inequahty 

the convergence in (fT^ also holds for the Hilbert-Schmidt norm, that is on the space L%. 
Using the symbol map and its inverse Op^ on T^r, we will convert operators on Lj^ 
into operators on L^(T^). We notice that for any A; G Z^, G for large enough N. The 
evolved density VuPk is in general not in X^r, but it is smooth since M is so. Any smooth 
density p is asymptotically equal to its projection on Tat, so that 

Vp G C°°(T2), WW"" o Op;v(p) - p||l2(t2) 0. 
Using this fact, Egorov's property can be recast into: 

VA; G I?, \\(yN{VM)pk - VuphWh-^iT^) 0. 

This means that the sequence of operators aN{T^M) semiclassically converges to Vm in the 
strong topology on Bi^L"^) (the bounded operators on L^(T^)). Now, a standard lemma in 
functional analysis ^ij Lemma 2.8] states that if a sequence An of bounded operators on 
some Banach space converge strongly to the operator A, then for any compact operator 
KAn converges to KA in operator norm. Since is compact, this implies that D^o(Jn{Vm) 
converges to V^^m in the operator norm. A simple comparison of the eigenvalues shows 
that ||criv(-De) — -Dg 11^(^2) — > as ^ oo, which achieves to prove the lemma. □ 

End of proof of the theorem: From lemma [H one applies standard methods to show 
that the spectrum oi a^iVe) converges to that ofVe, as was done in Section ITTl Namely, 
for any A 7^ and any small disk B{\,6) around it, the spectral projectors for aNiVe) 
and Ve in that disk satisfy ||n^i[^-| — nB(A,5)|| — > as A^ — > cxo. For small enough 5, the 
projector IiB{x,s) is independent of 6 and of finite rank, equal to the multiplicity of A in 
the spectrum ofVe- The above convergence implies that ^^^(xs) same rank for A^ 

large enough, and that the corresponding eigenspaces of Ve and a^iVe) are close to each 
other. Finally, the spectrum of is identical with that of a^iVe)- D 
This spectral stability was noticed numerically in 1^1] for the kicked rotator on the 
2-sphere. In their case, the coarse-graining consists in a sharp truncation of the Fourier 
expansion of the Husimi functions of the quantum densities, but the same arguments as 
above can be applied to show the spectral stability of the coarse-grained propagator in the 
semiclassical limit. 

4 On the non-commutativity of the semiclassical vs. small- 
noise limits 

In the previous section we have described the semiclassical limit of the quantum coarse- 
grained propagator for a fixed coarse-graining width e > 0. On the other hand, Section[21 
was dealing with the small-noise limit e ^ for the classical propagator Ve- 
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These two limits do not commute. Indeed, for fixed N the coarse-graining operator 
is a finite matrix, the eigenvalues of which converge to unity as e — (see Eq. (fT9|l ). 
Therefore, in this limit \\D^ — Idi-j^W ~^ 0- K{x) of compact support, we even have 
De = ^(^Ll; cis soon as the rescaled support {eN) Supp{K) has no intersection with Z^. This 
shows that for fixed and e decreasing to zero, is close to unitary on L^, uniformly 
with respect to the map M: 

^peLl s.t. IIpL, = 1, (l-||A-Id^^||)<||P,pL, <1. (22) 

On the opposite, there should exist a regime where N ^ oo (semiclassical) and simultane- 
ously e = e(A^) — * (vanishing noise) slowly enough, such that the eigenvalues ofVe{N) stay 
close to the eigenvalues ofVe{N), and therefore behave differently according to the classical 
properties of M. For an Anosov map, the "outer" eigenvalues (say, in some annulus Rr) 
will converge to the Ruelle resonances, while for an integrable map they will form dense 
strings touching the unit circle. Eq. (f22|l shows that a necessary condition for Ve to possess 
eigenvalues close to the origin (like V^) is that the coarse-graining operator itself has 
small eigenvalues. The smallest eigenvalues of correspond to the largest wave vectors 
in Z^, namely \k\ ^ N/2. From the explicit expression ^Ml, this implies the condition 

A^e(A^) » 1. (23) 

This condition is quite obvious: it means that the scale of coarse-graining e{N) must 
be larger than the "quantum scale" = 27ih, that is the distance between two nearby 
position states \qj)N- One may wonder if this condition is sufficient, or if e{N) should 
decrease slowlier to get the desired convergence. We have so far no definite answer to this 
question for a general map. 

In the next two subsections, we will compare the spectra of Ve and Ve for the various 
maps treated classically in sections 12.1.11 and 12.21 We know no nonlinear Anosov map for 
which Ruelle resonances can be computed analytically, so for this case we rely on numerical 
studies for the perturbed cat map [l]. 

We plot some numerically computed spectra, always using Gaussian noise and selecting 
two different ^-dependences for the noise width e(A^). We consider either a "slow decrease" 
e(A^) = A^~^/^, for which the convergence to classical eigenvalues is checked even for rel- 
atively small values of A^ (the largest value of A^ we considered is A^ = 40). To test the 
finer condition (f26|) . we also consider a "fast decrease" of the noise width e(A^) = ^"^^^ ^ 
the convergence to classical eigenvalues being then harder to verify numerically. 

We start by plotting on Fig. [T] (left) the spectrum of the quantum coarse-graining 
operator De. 

4.1 Quantum linear automorphisms 

In this section, we will only consider linear maps A satisfying the "checkerboard condition" 
given in Section [SIII21 necessary for their quantization on Hn- Then, the quantized linear 
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automorphism A satisfies the exact Egorov property |28| : 

Wk e I?, Afk/NA'^ = TAk/N VAOpNiPk) = OpNiVAPk) = OpN{pAk)- 

The quantum propagator therefore acts as a permutation on the quantum translations, hke 
the classical propagator on plane waves (cf. Eq. 0). In a first step, we treat any linear 
map, regardless of the nature of its dynamics. 

The main difference between the quantum and classical frameworks comes from the fact 
that Weyl quantization OpN is not one-to-one: OpN^Pk+Nk') oc OpN^Pk) for any k' G Z^. 
As a result, each orbit 0{k) = {A^k, t E Z} has to be taken modulo in the quantum 
case, yielding a "quantum orbit" ON^k) which is necessarily finite. Through a rescaling of 
j^, ON{k) is identified with a periodic orbit of the map A on the torus, situated on the 
"quantum lattice" (;^Z)^. The period T/v^ of this orbit is the smallest t > such that 
A^k = k mod Z^. To compute the spectrum of Va and Ve,A we need to analyze these 
quantum orbits. 




Figure 1: Spectrum of the quantum coarse-graining operator for the Gaussian noise (left); 
Spectrum of the coarse-grained 7r/2-rotation Ve,j (right). Parameters are N = 40, e = N^^/"^ 
(small circles). The large concentric circles are only shown for clarity. 

The quantum orbits form a partition of Z^, therefore a partition of the basis {Tk/N^ k G 
Z^} of L^. To each O^ik) corresponds the T/v,A:-dimensional subspace SpanON{k), in- 
variant through both Va and D^. The quantum propagator Va satisfies (we abbreviate 
T^^k with T) 

'PA.OpNipk) = OpNipATk) = OpN{pk+Nv) = (-1)'^ OpN{pk) (24) 

for some K G Z^, and ^ = k /\V + NV1V2. This implies that the eigenvalues of Va on 
SpanON{k) are the phases {exp(^(r + 7/2)), r = 0, . . . , T — 1}. By switching on the 
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noise, the equation (pljl is modified by inserting the action of on the successive modes 
TA^k/N- As a result, 



Vj'^y^fk/N - {-'^y d[^o(k) Tk/N, 



with 



T-l 

,(7V) def -TT ,(7V) 
i=0 



The spectrum of on SpanOp^{k) thus consists in T regularly spaced points on the circle 
of radius To estimate this radius, we take its logarithm 

^'ogMa,l = ^ E (25) 



fc'eciiv(fc) 



According to Eq. (fT9|l . this quantity is the average over the periodic orbit j^Oiy{k) of 

This function is smooth in ( except at possible 



the function fK{eN,() = log ^^(.jv.o) 
logarithmic singularities if 9k vanishes. In the Gaussian case G{x) = e-"l"l , /g has no 
singularities and admits for eA^ >> 1 the asymptotic behaviour /^(eiV, C) ~ ~7r(eA^C)^ iii 
the square {|Ci| < |, IC2I < |}- 

We have obtained an explicit relationship between the spectrum of V^^a and the struc- 
ture of periodic orbits of A on the quantum lattice. The latter drastically differs between 
chaotic and integrable automorphisms, which leads to qualitatively different quantum spec- 
tra. Below, we sketch the description of these quantum orbits, respectively for the elliptic, 
parabolic and hyperbolic maps. We use the notations and results of Sections 12.1.11 IT2^ 



4.1.1 Elliptic transformation J 

The quantization of the 7r/2-rotation J yields the finite Fourier transform J. The classical 
orbits 0{k) are of period 4. In general, the successive Pk are not congruent modA^, so 
that the quantum orbit O^ik) also has period 4. The only exception occurs for even, 
k = (y, y) (period 1) or k = (0, y) (period 2). Assuming that the coarse-graining kernel 
has the symmetry K{^) = K{JC,), the eigenvalues of j on a 4-dimensional SpanO^ik) 
are {i'(i[^'' , / = 0, . . . , 3}. Taking all quantum orbits into account, the spectrum forms 4 
strings as in the classical case (see Fig.[H right). For a Gaussian noise, Eq. TB^ shows that 
these eigenvalues are exponentially close to the corresponding classical eigenvalues. 



4.1.2 Parabolic linear shear S 

We can quantize on Hn the parabolic shear S described in Section 12.2.21 if the integer s 
is even. Quantizing the space of P^-invariants Vinv = Span{p(^ifi), I G Z} leads to the Vs- 
invariants Vinv,N = Span{T^j_ Q^, I G Z^}. Each T^^j^^Q^ is eigenstate of "P^ with eigenvalue 

'^1^0 z)' ■which is close to the classical eigenvalue K{e{OJ)) for eA^ >> 1; these eigenvalues 
then form a string connecting unity to the origin. 
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Figure 2: Spectrum of 'Pe,S2 foi' the linear shear S2 (circles) and values {(i^^j* q) i j^o (triangles). 
Parameters are = 20 (left), = 40 (right) and in both cases e = \og{N)/N . 

Now we study the spectrum of Ve on the orthogonal of the invariant space, V^^^^. 
For k = {ki,k2) with k2 ^ mod A^, the infinite orbit 0{k) = {k + (jsfc2,0), j e Z} 

projects modulo onto a finite orbit ON{k) whose period depends on the "step" g =^ 
gcd(Ar, sfca): Cjv(A;) = {(^i + jg mod A^, fca) | j = 0, . . . , ^ - 1}. For a fixed ^2 G Z*, the 
step g stays bounded in the limit of large A^, so that the sum (f25|) behaves as the integral 
Jq dt fx {eN, (t, ^)) ; for the Gaussian noise, the latter yields — 7re^(^ + fc^) + o(l). For a 
general value of A;2 G Z^r \ 0, the step g may be of order A^, in which case the sum (f25|l is 
not well-approximated by an integral; still, one can (for Gaussian noise) prove the uniform 
upper bound: 

ykeZ%\0, J-\og\df2^ \<-C{eN)' 

with C = 7rmin{l/s^, 1/16}. The eigenvalues of "P^ on V^^^^^r therefore situated on 
circles of radii < e"*-'*^^^)^ , so they uniformly converge to zero as eA^ 00 (we remind that 
the spectrum of Ves restricted to V^^^ reduces to {0}). 

In Fig 121 we show the spectra of Pe,52 fo'^ the linear shear 5*2 =( q ^ j and a "fast 

decreasing noise", together with the eigenvalues {d^^^^j q^}^^q. The 'non-invariant spectrum' 
converges slowly to the origin, while the 'invariant spectrum' becomes dense on the interval 
[0, 1]. Had we chosen a larger noise width, the non-invariant spectrum would be contained 
in a smaller disk for the same values of A^. 
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Figure 3: Spectrum of Ve^Ao for the cat map Aq. Parameters are = 30 (left), = 40 (right) 
and in both cases e = log(A^)/A^. = 30 corresponds to a "short quantum period" T(30) = 6, 
which is clearly visible on the shape of the spectrum. 

4.1.3 Hyperbolic transformations 

The space of invariant densities for a hyperbolic automorphism A reduces to Vinv = Cpo- 
Its orthogonal V^^^ = Lq is quantized into the space Lq ^ of traceless densities. We remind 
(cf. Section 1^.1. Ij) that the spectrum of Ve,A on this subspace reduces to the {0}, for any 
e > 0. 

The periodic orbits of a hyperbolic automorphism A were thoroughly studied in |37| : 
the authors classified the orbits according to the sublattice they belong to. Yet, 
they gave no equidistribution estimate on long periodic orbits. Our aim is to estimate 
the RHS of Eq. (f^ for all quantum orbits, at least for large enough N (we will restrict 
ourselves to the Gaussian coarse-graining). For any Anosov map, the long periodic orbits 
equidistribute in the statistical sense (averaging over all orbits of a given period) in the 
limit of long periods |3E1- For a certain class of hyperbolic automorphisms, semiclassical 
equidistribution of all quantum orbits was obtained for an infinite subsequence of values 
of [in]. Equidistribution morally implies that the sum behaves as the integral 
J^,dxfG{eN,x) ^ One can indeed show that for A^ in this subsequence, any 

eigenvalue A of Ve,A on the subspace satisfies 

|A|<exp|-^ + 0(e^Ar3/^)}. 

There also exist arbitrary large values of A^ for which the period T{N) of A modulo A^ 
is as small as + 0{1), where A > is the Lyapounov exponent of A All quantum 
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orbits ON^k) have periods dividing T{N). Starting from ko = (0, 1), the hnear dynamics 
shows that the point remains close to the origin (that is, at a distance << 1) along 
the unstable direction until the time ~ ^^j^, when it reaches the boundary of the square 
{|Ci| < 1/2, IC2I < 1/2}; it is then straight away 'captured' by the stable manifold, which 
brings is back to the origin during the remaining ~ steps. This orbit thus achieves 
an "optimally short" homoclinic excursion from the unstable origin, and is far from being 
equidistributed. The average of fc along this orbit is of order — ^j^^^, so the eigenvalues 

of Ve,A on SpanON{ko) have a radius ~ exp(— ^j^^^). These eigenvalues will therefore 
semiclassically converge to zero under the condition 

eN , , 

00. (26) 



We believe that these particular values of represent the "worst far as equidis- 

tribution of long orbits is concerned, and that condition (f^ suffices for the spectra of 
^e,yi|L2 to semiclassically converge to zero for the full sequence G N. In Fig. [HI we show 

the spectrum of V^^Aq for the quantization of the cat map = ( ^ with the fast 



decreasing noise width. 



Remark: For both the parabolic and the hyperbolic automorphisms, the spectrum of Ve 
on V{!^^ ^ reduces to {0} if the coarse-graining consists in a sharp truncation in Fourier 
space, and grows as ~ cN with c a finite but small constant (depending on the 
classical map). Any nontrivial quantum orbit ON{ko) then contains an element k s.t. 
\k\ > e~\ such that the corresponding mode Tk/N is killed by D^. 

4.1.4 Quantum translations 

As we mentioned in Section l3.1.2[ any classical translation Ty with f G can be quantized 
on TYat by the quantum translation T^(jv), where the "closest" point to v on the 

quantum lattice. This quantization was shown [22] to satisfy Egorov's property (fT4ll . 

From Eq. (fH^ . the quantum propagator V^iN) = ad(T^(jv)) admits any quantum trans- 
lation Tk/N = OpN{pk) as an eigenstate, with eigenvalue e"^''^^*'^'^'^. All eigenvalues are 
A^-th roots or unity, and are at least A^-degenerate; in case the classical translation is 
ergodic, these degeneracies contrast with the nondegenerate (but dense) spectrum of Vy. 
On the opposite, for a rational translation (the classical eigenvalues form a finite set), the 
quantum eigenvalues may take values in all A^-th roots of unity, in case A^ and v'^^'^ are 
coprime. The spectra of Vv and P^(jv) may thus be qualitatively very different. 

This difference disappears when one switches on the noise. Each Tk/N is also an 
eigenstate of the coarse-grained propagator V^^v = D^Vy(N) , associated to the eigenvalue 
g2i7rkAv(f^)|^(N) ^j^^ scmiclassical limit, the deviation from the corresponding classical 
eigenvalue e^^'^'*^^K(ek) (cf. Section 12.2.111 depends on both the wavevector k and the 
difference v — v^^\ 
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Figure 4: Spectrum of Ve,v for the translations v = (0, 1/3) (left: N = 30, center: N = 37) and 
V = (l/-v/2, l/\/5) (right: N = 37). The noise strength is e = A^~^/^. For the rational translation, 
the eigenvalues are either exactly on the classical axes (for N multiple of 3), or semiclassically 
converge to it. For the irrational one, the eigenvalues become dense in the full disk. 

To give an example, the rational translation vector v = (0, |) leads to 3 eigenvalues 
g2i7ri/3 Jqj, ^Yie classical propagator P^, and 3 (infinite) strings for the spectrum of "P^^e. 
Quantum-mechanically, if is a multiple of 3 one takes and the eigenvalues of 

Vy^e are exponentially close to the classical ones (Fig. lU left). In the case = 3n + 1, 
the quantum translation vector v^^ill be v^'^^ = (0, 3^^), so that the classical and quantum 
noiseless eigenvalues for the mode pk deviate of an angle 2Trk A {v — v^^^) = ^f^- this 
deviation can be as large as tc/S for wavevectors \ki\ ^ N/2. After switching on the noise, 
the classical and quantum eigenvalues for k G can deviate by at most 0{-^), the 
maximal deviations occurring for wavevectors |A;i| ~ e^^, k2 = 0{1) (Fig. |31 center). 

4.2 Examples of quantized nonlinear maps 

In this section we shortly review the spectrum of quantized coarse-grained propagators 
for three nonlinear maps: we first treat the nonlinear parabolic shear of section 12.2.31 and 
the stroboscopic map of the Harper Hamiltonian (Section 12.2.4^ . for which we have some 
analytic handle. We then consider a perturbation of the cat map Aq and compare the 
numerically computed eigenvalues of Ve to the spectrum of obtained in [T|. 

• the nonlinear shear described in Section 1^.2.31 is quantized by taking the product of 
the quantum linear shear S (see Section I4.1.2jl with the matrix e~'^/'^. Since F{p) 
is a function of the impulsion only, its Weyl quantization F acts diagonally on the 
impulsion basis As a result, the perturbation ad(e~^^/'*^) acts trivially on any 

projector \pj){pj\, and therefore on any translation T(o,m)/Ar- As a result, the spectrum 
of the noisy nonlinear propagator on Vinv,N forms the same string as for the linear 
shear. 

The action of V on V^^^ ^ can be studied as in the classical case: the propagator 
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acts inside each subspace Vn^N = Span{T(n,m)/N, ^ £ '^n}, as a unitary N x N 
Toeplitz matrix V^"'^ (instead of a simple permutation). In Appendix IA.2I we study 
the non-unitary spectrum of Pe"^ when taking for K{^) a sharp cutoff: Pi"'' is then 
the truncation of P*^"^ to the subspace {|m| < e~^}. We compare this truncated prop- 
agator with the corresponding classical one, and show that if eN >> 1, both spectra 
belong to the same union of 1-dimensional strings contained in a fixed 'small' disk 
around the origin, the size of which depends on the strength of the perturbation F'. 
The spectrum of V for the quantized nonlinear shear e~'^/'^S2 with F = cos(27rp) 
is shown in Fig. El (left) for = 40 and e = log(A^)/A^, together with the values 

t".,(i,o)/i=o- 

The quantization of the stroboscopic map ip]^ for the Harper fiow, i.e. the unitary 
matrix e~^^^^, leads to the propagator Vh leaving invariant any density of the type 
H"^. Under the same conditions for the coarse-graining kernel K{^) as in the classical 
case (see Section l2.2.4|) . VH,e may admit for eigenstates H, {H^ — Id) and (ff^ — 
7+2cos(27r/Af) ^-^ , and the associated eigenvalues {(i[^jQ)}^=i converge to the classical 
eigenvalues {k{je)}^^^ if eA^ — > oo. The spectrum of VH,e for A^ = 40 and e = N'^^^ 
is shown in Fig. [HI (right), together with the values {(i^^jo)}i=o- For the Gaussian 
kernel used there, only H is an eigenstate of Z)^, so the subsequent eigenvalues are 
hybridized by the coarse-graining. A "string" of eigenvalues along the positive real 
axis is clearly visible on the plot. 

We consider the quantization of the perturbed cat map studied in [T], that is the 
map : ^ (3g+2p+.(coS)-cos{4.g))) • We selected the perturbation k = 0.5/27r 
in order to compare the quantum spectrum with the classical resonance spectrum 
described in [H Fig. 2]. In Figure El we plot the spectrum of for A^ = 40 together 
with the 7 'outer' resonances whose values are given in [H Table I]. In the 'large noise' 
regime (left), the largest quantum eigenvalues are indeed close to these resonances, 
and the rest of the spectrum inside a smaller disk. The small-noise regime (right) does 
not uncover the classical resonances, however the non-invariant spectrum is already 
contained in a relatively small disk around the origin, of radius comparable with the 
largest resonance. 



Concluding remarks 

In spite of the relatively low values of A^ used in the numerical plots, one can already 
clearly distinguish two different behaviours of the quantum spectra, depending on the 
classical motion, especially for the slow noise decrease e = A^~^/^. For the integrable maps, 
the largest nontrivial eigenvalue of the noisy quantum propagator is at a distance from 
unity, and this eigenvalue is the first one of a string of eigenvalues connecting unity with a 
neighbourhood of zero. On the opposite, for the two chaotic maps we studied (linear and 
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Figure 5: Spectrum of Ve for the nonlinear shear e~'^/'*S2 (left, A?^ = 40, e = \og{N)/N) and the 
Harper map c~^^/^ (right, N = 40, e = N~^/'^). In both cases we also plotted some values d^f^^Q^ 
(black triangles). 
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Figure 6: Spectrum of Pe for the quantum perturbed cat map A^i (empty circles) together with 
the largest 7 classical resonances (black diamonds). Left: N = AO, e = iV-V2. Right: = 40, 
e = log{N)/N. 
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perturbed cat maps), the spectrum already shows a finite gap for these values of A^, even 
in the fast-decreasing noise regime e = log(A^)/A^; however, for the relatively low values of 
N shown in the plots, the gap is governed by the classical resonances only in the regime 
of larger noise e = N'^^"^. I conjecture that the quantum spectrum also converges to the 
classical spectrum (in any annulus Rr) in the regime e = log(iV)/iV, this being visible only 
for much higher values of A^. 

General conclusion 

In this paper we have shown the connection between the spectra of coarse-grained quantum 
and classical propagators for maps on the 2-dimensional torus. I claim that Theorem [T] 
can be straightforwardly extended to maps on any compact phase space, like the 2-sphere 
5^ used in (Hll EE]- Using our knowledge of the spectrum of the classical noisy map, one 
infers the presence or the absence of a finite gap between the (trivial) eigenvalue unity on 
the one hand, the rest of the spectrum on the other hand, taking both the semiclassical 
{h — >■ 0) and small noise (e — ^ 0) limits in a well-defined way. The presence of this gap in 
the classical noisy propagator is related with the ergodic properties of the map. On one 
extreme, corresponding to integrable maps with infinitely many invariant densities, there 
is no gap in the small-noise limit, and the spectrum contains a dense 'string' of eigenvalues 
connecting unity with the interior of the disk: the long-time evolution is therefore governed 
by the eigenmodes associated with these large eigenvalues, and is typically of diffusive type. 
On the opposite extreme, for an Anosov map ("strongly chaotic"), the spectrum exhibits 
a finite gap, responsible for the (classical) exponential mixing. The case of maps with less 
chaotic behaviour is far from being settled. The example of irrational translations shows 
that ergodicity does not imply the presence of a gap. Exponential mixing was recently 
proven for some non-uniformly or partially hyperbolic maps (see |0| for a review of recent 
results), which should (?) imply a gap in the spectrum of V^- Maps with subexponential 
mixing would also deserve to be studied from this point of view. The more general case 
of maps with mixed dynamics (where the phase space can be divided between "regular 
islands" and a "chaotic sea") was mostly investigated numerically [211 EI]; one found the 
presence of both "dense strings" of eigenvalues and isolated resonances in the spectrum of 
Pg. This (physically relevant) type of dynamics certainly deserves further study. 

Like in the classical framework, the presence of a gap in the spectrum of the quan- 
tum noisy propagator allows to describe the long-time dynamics of densities induced by 
this propagator: p(t) = ('Pe)*p(0). This dynamics is a possible model for a dissipative 
perturbation of the unitary evolution p M^pM~^, the dissipation being induced by the 
interaction with the environment [7J. One measure of the decoherence occurring through 
this evolution is the "purity" tr(/)(t))^. For an Anosov map M, the long-time evolution 
of this quantity will be governed by the largest nontrivial eigenvalues of V^, which are 
semiclassically close to the classical resonances. This very evolution was recently studied 
for perturbed quantum cat maps, providing an estimate of the largest eigenvalues of for 
large values of N [21]. Another way to characterize the time evolution of Ve is through the 



30 



dissipation time tdiss = infjt G < e~^}. In the case of linear automorphisms 

on the (i-dimensional torus, the small-noise behaviour of tdiss has been shown to qualita- 
tively depend on the dynamics jlHj. The same conclusion should apply to non-linear maps 
as well. 

Our noise operator consists in an average over random translations in phase space. 
To diagonalize this operator, we have used the fact that the classical (resp. quantum) 
propagators for translations Vv (resp. Vv) form a commutative algebra, and admit as 
eigenfunctions the Fourier modes. We obtained explicit expressions for the spectrum of 
Ve^M for the class of linear maps M, because these maps acted simply on these Fourier 
modes. However, as noticed in Section one may be interested in a more general 

type of diffusion operator, like the one defined in Eq. (pT|l . The spectrum of Dj^^^ o Vm 
for a given map M should be qualitatively independent on the family of Hamiltonians 
{Hj}, as long as the operator De ' has the same characteristics as the coarse-graining D^, 
namely it leaves almost invariant some "soft modes", while killing fast-oscillating modes. 
In particular, for M an Anosov map, the outer spectrum of ^ o Vm should converge to 
the set of resonances {Aj}, and the rest be contained in a smaller disk, in the limit e — 0. 
For an integrable map, the eigenvalues will not pointwise converge to the ones of DJ^m-, 
but they should also accumulate along a 'string' touching unity. 

Our choice for the superoperator was inspired by its classical analogue, as well 
as its relevance to modelize the "quantum noise". However, in some quantum systems 
experiencing (weak) interaction with the environment, the effective "noise" or "decoherence" 
superoperator may have no obvious classical analogue; this seems to be the case for instance 
in Nuclear Magnetic Resonance experiments aiming at realizing a "quantum computer" 
[T5] . The quantum Hilbert space is composed of a sequence of n spins, so it is isomorphic 
with the torus Hilbert space TYat for N = (each spin corresponds to a binary digit of 
the position qj). The decoherence operator acting on such a system (for a small number 
of spins) has been experimentally measured, it resembles the product of operators acting 
independently on the individual spins, and therefore does not enter in the family of diffusion 
operators described in the present article; in particular, such a decoherence operator seems 
to have no obvious classical counterpart. Conjugating this type of decoherence operator 
with a unitary evolution (typically, the finite Fourier transform J) may lead to a drastically 
different spectrum from the one shown in Fig. [T] (right) , not excluding the presence of a 
gap. 
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A Appendix 



A.l Classical nonlinear shear: spectrum of truncated Toeplitz ma- 
trices 

This appendix describes the spectrum of the noisy classical propagator for a nonlinear 
shear fSection l2.2.3p . In each invariant subspace Vn = {e^'™''p(p)} (n G Z*), the Perron- 
Frobenius operator V acts as a multiplication of the function p{p) on the circle 3 t = 

a„(e2-P) Q-2iMsp+F'(p)) ^ J2 &„(m)e2i™P. (27) 

Notice that a„(t) = (ai(t))". In the Fourier basis, this multiplication has the form of an 
infinite Laurent matrix -L(a„), with entries L(a„)jj = a„(z — j). From standard results, 
the spectrum of V\v„ is absolutely continuous, with support the range of the function a„(t) 
(here, the unit circle). 

We consider the coarse-graining consisting in a sharp truncation in Fourier space: 
K{C) = 0(^i)0('C2)- The noisy propagator Ve then acts on Vn as the matrix L(a„) trun- 
cated to the block {\m\, \m'\ < e^^}, or equivalently as {1 < m,m' < 2e^^ + l}. We assume 
that the smooth function a„(t) on can be continued to an analytic (or meromorphic) 
function on C* (for the linear shear, one has an{z) = 2;"'*"). Then, for any r > and any 
e > 0, the spectrum of this truncated matrix is contained in the convex hull of the curve 
an(rT^) Prop. 2.17 and Section 5.8]. For the linear shear, these curves are the circles 
centered at the origin, so the spectrum reduces to {0}, as expected. If the perturbation F' 
is small, these curves are deformations of circles, some of which remain in a small neigh- 
bourhood of the origin. As a concrete example, the perturbation F{p) = ^ cos(27rp) leads 
to the function a^i{z) = z'^ exp {na^l/ z — z)}. The curves a_i(rT^) satisfy 

Vr > 0, = r =^ |a_i(2;)| < exp(7ra|l/r — t|). 

We assume that s > 0. For a small perturbation (a << 1), the function on the RHS has 
a minimum of value amm(a) ^ (^) ■ This implies that for any n, the spectrum of Ve\Vn 
is contained in the disk of radius (amj„(a))'"' for any e > (we have used the symmetry 

Cl—n C-n)- 

For a more general perturbation 2nF'(t) = J2m>i fmt"^ + c.c. one can use the bound 

27r|3F'(rt)| < sup(5'(r), g{l/r)) with g{r) =^ J2m>i l/mk™- If "small", then the 

function exp{g{l/r)) has a minimum a^m << 1 on M.'^. This implies that the spectrum 
of Ve\v„ is contained, for any e, in the disk around the origin of radius a[^|„. 

A. 2 Quantum nonlinear shear 

The results of the previous appendix can be easily adapted to the quantized nonlinear 
shear. The latter is defined as the composition of the quantum linear shear S with the 
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Floquet operator corresponding to the Hamiltonian F (see Section l2.2.3j) : the quantum 
map reads Se~'^^'^^'^ . Its action is diagonal in the impulsion basis {\pj)N} ofHN {pj = j^, 
with j G Ziv): 

^e-FA|pj)N = e-2^-^(tP.'+^(Pj)) |pj)N. 

Since s is even, this expression is well-defined for any j G Z^r. As we now show, it allows 
to express the action of the propagator V on the quantum translations T^/n- Taking 
k = {m,n) G Z^, this translation can be decomposed as 

V acts inside each subspace Vn^N = Span{T(^rn,n)/N , ^ ^ '^n} = Span{\pj+n){Pj\, j ^ '^n}- 
The propagator V indeed multiplies each \pj+n){Pj\ by the phase A^^at (e^''^^Pj^2N)^ ^ with 
the function An^^ on defined as: 

An,N{t) uniformly converges to the function an(t) of Eq. (f27|l in the limit ^ oo. For our 
example F{p) = ^cos(27rp), it takes the concise form An,N{e'^'''''^) = e-2i™{«P-°".N sm(27rp)) 
with an,N = a£ sin(7rn/Ar) = «(1 + OUn/Nf)). 

From there, we can express the action of V on the quantum translations: 

V T(m,n) = f(m-l.n) y4„ Ar(/) = \^ T (i,n) A'^ J^^TTI — I), 

N N N ' 

where A'^]^{m) = X]jyGz(~-'-)"''^«.A^("^ + uN). Therefore, V acts on Vn,N through the 
N X N Toeplitz matrix with coefficients y4^^(m). If we truncate this matrix to the size 
{1^1, l^'l < with e^^ << A^, we only take into account the coefficients with \m\ < 2e^^. 
Since An^^ is a smooth function, each of these coefficients is the sum of a "dominant" term 
An,N{'m') and a "remainder". For our example, the classical coefficients are given by Bessel 
functions: a„(m) = Jm+sn(2vr?T,a), and the same for A„ jv(m), replacing a by an,N- These 
coefficients therefore satisfy |74„ Ar(m)| < C' ^']°t+ns|! ' ' ^^^^ remainder is uniformly 
bounded above by (7^)^^^ for large A^. If we call V^^y ^ the truncated Toeplitz matrix 
with coefficients A„^iv(m), we get the estimate: 

This estimate implies that the spectra of both matrices cannot be at a distance larger than 
0{{eN)~''^/^). The matrix Vj^y ^ may be analyzed along the same lines as Ve\v„ in the 
previous Appendix. Its spectrum is contained for any e > and any r > in the convex 
hull of y4„ Ar(rT^), which converges to the convex hull of a„(rT^) when A^ ^ 00. 
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